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Abstract 

Mathematical model that allows for direct tracking of the homoepitaxial crystal 
growth out of the window etched in the solid, pre-deposited layer on the substrate is 
described. The growth is governed by the normal (to the crystal-vapor interface) flux 
from the vapor phase and by the interface diffusion. The model accounts for possibly 
inhomogeneous energy of the mask surface and for strong anisotropies of crystal-vapor 
interfacial energy and kinetic mobility. Results demonstrate that the motion of the 
crystal-mask contact line slows down abruptly as radius of curvature of the mask edge 
approaches zero. Numerical procedure is suggested to overcome difficulties associ- 
ated with ill-posedness of the evolution problem for the interface with strong energy 
anisotropy. 

Keywords: Thin films, epitaxy, MOCVD, surface diffusion, interface dynamics, 
contact lines, rough surfaces, wetting, regularization of ill-posed evolution problems. 
PACS codes: 68.10.Cr, 68.35.Ja, 68.45.Gd, 68.55.-a, 81.10.Aj. 

1 Introduction 

The Epitaxial Lateral Overgrowth (ELO) and Selective Area Epitaxy (SAG) are commonly 
used to grow micro-scale semiconductor crystals and thin films. The micrometer-scale selec- 
tive growth is well-achievable by chemical vapor deposition or liquid phase epitaxy [l]-[4]; 
the ELO and SAG can be used also to grow nanostructures by the molecular beam epitaxy 
[5]-[7]. 

Recently, the mathematical models that allow to numerically study the ELO and SAG 
from vapor [8] - [10] were introduced. These continuum, geometric crystal growth mod- 
els are formulated as free-boundary problems in two dimensions (that is, normal to the 
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substrate/mask and normal to the axis of long stripe openings etched in the mask) and, 
therefore, they are capable of explicitly accounting for the mask topography and the inter- 
face anisotropy. The crystal growth is by the normal flux of particles from the vapor phase 
to the interface; the surface diffusion redistributes the material along the interface during the 
growth. Evaporation-recondensation from both the mask surface and crystal- vapor interface, 
and the diffusion of adatoms along the mask surface to the tri-junction crystal-mask-vapor 
are accounted for by the models. The details of the formulation, boundary conditions and 
numerical method can be found in [8, 9]. In this paper we only address issues related to the 
mask surface modeling and the crystal interaction with the edge of the mask (mask corner 
point). Note that the crystal growth out of the trench occurs relatively early in the ELO 
process; however, our previous modeling suggests [8, 9] that the perturbation introduced to 
the crystal-vapor interface by the edge of the mask exists for a long time after the actual 
interaction of the advancing edge of a crystal (contact line) and the mask edge took place. 
Thus it seems important to understand the details of the contact line motion over the edge. 

In [8, 9] the treatment of the contact line motion relies on the assumptions of (a) the 
thermodynamic equilibrium at the tri-junction (indeed, SAG and ELO are equilibrium pro- 
cesses), (b) the mask surface as the rectangular, macroscopic step, and (c) the mask as 
the energetically homogeneous surface: the mask-vapor surface energy / ~f mv = const. The 
assumptions (b) and (c) limit the ability of any model to describe the overgrowth of the 
crystal line from the vertical onto the horizontal part of the mask. Indeed, (a) implies that 
the equilibrium contact angle is given by the equation 

Imv -1cm -7«;COS0+ SHI = 0, (1.1) 

where 7 cm and 7 C1 , are the energies of the crystal-mask and crystal-vapor interfaces, respec- 
tively (see Appendix A). Then, (c) together with the natural assumptions of constant (at 
the tri-junction) 7 cm and 7 C „ imply that is a constant value at the mask. However, due to 
(b) the mask edge is a singular point in the sense that the contact angle is not well defined 
there. In [9], this situation was dealt with by instant readjustment of the interface — x-axis 
angle (angle of the tangent, <p x ) to a new value once the contact point reaches the mask edge 
to preserve the equilibrium constant angle 0. Such instant readjustment causes the short 
unphysical retraction of the contact point into the vapor phase and influences the shape 
of the interface x . In this paper we make an attempt to model the overgrowth event self- 
consistently by approximating the mask edge by a circular arc and restricting the contact 
point to movement along this arc. We numerically study the limit of zero arc radius. Since 
in applications the mask is often a crystalline solid itself, we allow for the anisotropy of i mv 
and thus, by virtue of (1.1), for the variable along the mask. (Quartz or tungsten masks 
are common [2]. Sometimes a mask is not deposited at all, but the long stripes are etched 

1 Angle of the tangent is necessary for the determination of the location of the contact point on the mask 
by coordinate's extrapolation from marker particles adjacent to the contact point, through 

dy/ds 

where s is the arc length along the curve (crystal- vapor interface). 
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in the crystalline GaAs or sapphire substrate itself [1, 11, 12]. The long crystalline mesas 
left on a surface of a substrate in this latter case are more advantageous for the epitaxy of 
some semiconductor materials than the geometrically equivalent mesas of the foreign (to the 
substrate) mask material.) 

Apart from the interest to this problem that stems from the ELO technology, there is a 
broad class of phenomena in which the material-defect interaction is a key component. For 
instance, when liquid wets the microscopically rough surface the contact line gets pinned on 
surface asperities and thus becomes unstable. The latter problem, both for a single asperity 
and for the random ensemble of asperities has been studied intensively from a physico- 
mechanical standpoint; see, for instance, recent works [13, 14] and the review article by de 
Gennes [15] . The somewhat different situation where the surface defect is the macroscopic, 
organized structure was considered, for example, by Oliver et al. in [16]. In that work, the 
defect is the edge of a solid disk that supports the liquid drop on its top. The pinning effect of 
the edge has been examined theoretically and Gibbs inequality condition for the equilibrium 
of a drop bound by the edge has been confirmed experimentally. Also, the various conditions 
for drop stability at the edge have been examined. Authors used a circular sapphire disk 
with the 90° edge and the aluminum disks with the edges subtending a range of angles. In 
[17], the finite element stability analysis of an inclined pendant drop was performed, as well 
as the experiment. The latter showed that both the contact line and the contact angle can 
adjust around the capillary until the Young-Laplace equation that balances the capillary 
pressure with pressure forces due to gravity [15] is satisfied. In [18], the asymptotic solution 
to the model of the contact line pinning on a single, macroscopic spherical defect was de- 
rived. However, our bibliography search for the study in which the direct numerical tracking 
of the contact line motion over the defect is performed was unsuccessful. Obviously, the ELO 
problem for the homoepitaxial crystal is simpler than any similar problem involving a liquid, 
since there is at least no need to solve the Young-Laplace equation. Heteroepitaxial and 
thermal stresses, if present in a solid-on-solid system, may significantly complicate the anal- 
ysis, but they usually are incapable of modifying the equilibrium contact angle determined 
by the interfacial energies alone [19] (however, as is well known, they may alter the dynamics 
of the free surface of a solid film) . The latter in the case of the SAG was demonstrated in 
[20] , where the morphological evolution of a heteroepitaxial thin film growing on a patterned 
substrate was studied by the phase-field method. 

2 Formulation 

The Figure 1 shows the problem geometry. As was already noted, the length of the open, 
unmasked stripes is assumed much larger than their widths since large aspect ratios are 
common in the applications; thus, the 2D crystal growth model can be used. 
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Figure 1: A sketch of the mathematical situation for the crystal growth on the stripe- 
patterned, masked substrate. The crystal-vapor interface is defined parametrically as y — 
y(s, t), x = x(s, t), < s < S(t), where s is the arc length along the curve and S is the total 
arc length of the curve. Due to periodical arrangement of the stripes, the crystal growth is 
studied on a partial cross-section which is a line segment extending from the center line of 
one mask surface at x = — £ to the center of the adjacent stripe at x — L. The interface 
is sketched such that the crystal overgrowth on the mask is shown. The mask is assumed 
rectangular, with the rounded edge (solid curve; since the lateral length scale in the picture 
is much larger than the vertical one, the edge in the form of the circular arc does not look 
circular), a is the angle that the outward unit normal to the interface, q c , makes with 
the horizontal axis. is the contact angle that the interface makes with the mask surface, 
measured from within the crystal. 9 is the angle that the unit normal to the mask surface, 
q m , makes with the horizontal axis. (x*,y*) is the contact point. 
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2.1 Equations of crystal growth 



Based on the classical Mullins' approach [21], the normal velocity of the interface is formu- 
lated as a function of the local curvature and its derivatives. In this work, we use 



V a = D 



ds 2 



Jcv{u) + 



d 2 % v {a) 
da 2 



K 



+ JM(a), 



(2.1) 



which is the simplified expression from [9] . We neglect the evaporation- recondensation term 
as being unimportant for this study and assume the isotropic interfacial diffusivity. The 
equation (2.1) is in the non-dimensional form. D has the meaning of the constant effec- 
tive interfacial diffusivity (see Table ?? in the Appendix B), K is the curvature and J the 
composite parameter that involves, in particular, the chemical potential of the vapor phase 
and the mean value of the interface mobility, Mq (the kinetic coefficient), ^(a) and M(a) 
are the anisotropy factors of the interface energy and mobility, respectively. These factors 
are specified in Section 2.4. The marker particles on the interface are numerically advanced 
using the parametric evolution equations for the Cartesian coordinates of particles, viz. 



dx 
dt 



V„ 



-v„ 



dy_ 

ds ' 

dx 



(2.2) 



dy_ 

dt Xq ds' 

These equations are fourth-order parabolic when V q is given by (2.1). 



2.2 Mask Energy and Shape 

To derive the expression for 7 mt) which accounts for the additional edge energy, we use the 
method of Xin & Wong [22, 23], which they employed in [24] for studies of the grain boundary 
(GB) grooving in polycrystalline thin films. In contrast to the case of dynamically evolving, 
faceted interfaces such as GB grooves, the application to the static crystalline mask surface 
is straightforward since the method requires an a priori knowledge of facets orientations and 
lengths. Details of the method can be found in [22, 23]. For the rectangular crystal in the 
thermodynamic equilibrium that represents the mask on a substrate (half of that crystal 
is shown in Fig. 1,2), the result is the following nondimensional expression (units for the 
surface energy are chosen ^ m h m /Q, where Q is the atomic volume and \i m is the equilibrium 
chemical potential of the mask surface): 



e-r + 2(r-e)F[e-- 



cos# + 2(l -r) 



sin#, 0<9<n. 

(2.3) 

In (2.3), r = Rj H is the ratio of the nondimensional circular arc radius to the nondimensional 
step height, t — d/H, where d is the nondimensional half-width of the mask, and F(u) is 
Heaviside function (all lengths are nondimensionalized by the half-width of the stripe, L). 
For our modeling of the ELO, we are only interested in values of 6 in [0, 7r/2]. The function 
7„™(6 ) ) is plotted in Fig. 3 for I = 1 and decreasing values of r. Notice that ^ mv is at 
maximum at 6 = ir/4, since the half-mask is taken square with the rounded edge. The 
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sharper is the edge, the larger is the maximum value. For the 90° edge r is exactly zero, 6 
takes on values and tc/2 only, and (2.3) gives 7 mu (0) = r y mv (n/2) = I. 




Figure 2: Half of the mask with the rounded edge. 
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Figure 3: The nondimensional energy of the mask surface as function of the angle of the unit 
normal (to that surface), for the different ratios of the radius of the circular arc to the mask 
height. The points (0, 1) and (vr/2, 1) correspond to the vertical and horizontal parts of the 
mask, respectively; other points correspond to the circular arc that connects the vertical and 
horizontal parts. 



To parametrize the mask surface, we take angle 6 G —Oi,02 + tt/2 as the independent 
variable (see Fig. 2). From geometry, 



Then, 



n H-R (\ \ 

V\ = arctan — - — = arctan 1 , 

R \r J 

~ d-R (d \ 

02 = arctan — - — = arctan I — — 1 I . 

R \ R I 



o, if - #i < e < o, 

8, if < < tt/2, 

tt/2, if tt/2 < 6 < tt/2 + 2 . 



(2.4) 



(2.5) 
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Also, 



o, if - e 1 < e < o, 



+ cos#) , if < < tt/2, 
-R (l + tanfl) , if vr/2 < < tt/2 + 2 , 

' H - R (l + tanfl) , if - ft < < 0, 
y m = < if - i2 (l - sin fl) , if < 6> < tt/2, 
k if, if n/2<9< tt/2 + # 2 



(2.6) 



(2.7) 



are the nondimensional parametric functions that describe the mask surface. From (2.5)- 
(2.7), the rounded portion of the mask is described by 




ff(-l + cos0), 

H - R (1 - sin 9) , < 9 = < tt/2. 



(2.8) 



Taking /i m h m /fl for the unit of 7 cm , and 70 (the mean value, see (2.16) below) for the 
unit of 7c„, the equation (1.1) takes the form 



Eijmv - 1cm) - %v(a*) cos0 + "^(«*) sin (f> = 0, 



(2.9) 



where E = fi m h m / (^7o), Imv is given by (2.3), 7 cm now is the nondimensional energy of 
the crystal- mask interface (assumed constant and isotropic), ^ cv is the anisotropy factor of 
the interfacial energy, and a* is value of the angle of the normal to the interface at the 
tri-junction. Notice that 

a* = (f) + d. (2.10) 

(2.9) is solved once in the beginning of the computation and yields values of the imposed 
contact angle at every 0i e —0i, 02 + ^/2 ,i = l---Ng, where Ng is the number of equidistant 
nodes in the interval. Ng is taken in the range 500. ..3000, the larger the smaller is the radius 
of the circular arc (since for the accurate marching of the contact point along the arc it is 
important to have sufficient number of nodes there). 



2.3 Updating the contact point 

In this section, the procedure which updates the position of the contact point on the mask is 
described. The update is needed after all marker particles on the interface but the contact 
point are moved one time step forward using the evolution equations (2.2). 

Let x* and are the "old" coordinates of the contact point. Using x*,y*, the coordinates 
of the two nearest nodes on the mask: x m {i) , y m (i) , x m (i + l),y m (i + 1), and 0^+1, <fii,i+i, we 
first compute values of 9 and at (a;*,?/*) by linear interpolation. Using the latter values, 
the angle of the tangent to the interface at (x*, y*) is found from 

0^ = + ^-^/2. (2.11) 
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Let (xi,yi) and (x 2 ,y2) are the coordinates of the first two marker particles on the interface 
closest to the contact point. In case the contact point is on the vertical part of the mask, 
the updated location is (0,2/*), where 

V* = \ (%i -V2 + ( x 2 ~ 4xi)tan0 x ) . (2.12) 

This follows from the second-order, one-sided finite difference approximation of the contact 
angle condition (1.2). In case the contact point is on the horizontal part of the mask, the 
updated location is (x*, H), where 

t. = l( iXl - Xl + 3H - 4 * + » ). (2.13) 

3 V tan^ J 

In case the contact point is on the circular arc, we solve the quadratic equation (ref. (2.8)) 

I (4yi - y 2 + (3x* - 4xi + x 2 ) tan^) =H-R(l- yfl - (1 + x*/R) 2 ^ (2.14) 

and find two ^-coordinates x^\x^ and corresponding y-coordinates y* 1 ^, yi 2 \ Then, if both 
x£\x*P > —R (that is, if both possible updated locations are on the arc), we take the pair 
of coordinates closest to (x*, y*) for the updated location. If only one of (xi 1 , y* 1 ^), (x* 2 \ y* 2 ^) 
is on the arc and the other is not, we choose the pair that is on the arc (closer to the "old" 
location). 



2.4 Anisotropy 

The anisotropy factor of the kinetic mobility is chosen as in [25], where the following form 
is suggested and the simulations of faceted growth are performed within the framework of 
phase field bulk solidification model: 

M(a) = 1 -a + 2at&nh- — -— ; M = M M(a). (2.15) 

V ; |tan2(a + /? M )| 

In (2.15), a and \ are constant parameters that control the width and flatness of facets, 
respectively; (3m is the phase angle. Facets are formed at regions on the interface where 
a = ir/4 — (3 M + nir/2, n = 0, 1, .... It was shown in [25] that the results of computations with 
the four-fold anisotropy (2.15) are in good agreement with the so-called kinetic Wulff shapes 
(see references in [25]; specifically for the selective area epitaxy the kinetic Wulff shapes were 
constructed and compared to the experimental shapes in [26]). For the simulations discussed 
in the next section we chose a = 0.95, x — 2, Pm = 0. 

The anisotropy factor of the crystal-vapor interfacial energy has standard, four-fold sym- 
metric form 

% v (a) = 1 + e 7 cos [4 (a + /3 7 )]; ^ cv = 7o7«>(a), (2.16) 

where the constant e 7 > determines the degree of anisotropy, and (3 7 is the phase angle. 
In this work we allow for strong anisotropy of 7 C „, e.g. e 7 > 1/15. In this case, the interface 
stiffness G = ^ C v + d 2 ^ cv / da 2 is negative in certain intervals of values of a, and that manifests 
in the appearance of corners on the equilibrium Wulff shape of a crystal [27]; in the dynamical 
case, the evolution equations become backward-parabolic and unstable where G < 0, see for 
example [28] - [30]. 
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2.4.1 Regularization 

To track interface evolution with the strongly anisotropic interfacial energy, the problem must 
be regularized in order to penalize spatial oscillations as well as the tendency to form corners. 
The regularization we use in this work is based on the addition of curvature dependence to 
7 C „. Regularization by curvature was first proposed by Herring [27] for the equilibrium 
shape; see [31] for the modern mathematical analysis of this problem. Apart from being 
useful as a modeling tool, the regularization by curvature has physical ground, namely the 
interaction of steps on a crystal surface in the vicinity of a corner [29]. Unlike [29, 30] and 
other works where the 5 k 2 , 5 = const, (k is the dimensional curvature) term is added to the 
surface energy on the problem formulation stage, here we regard the regularization only as a 
numerical technique intended to keep the computation alive, and regularize by adding same 
term to 7 C1 , given by (2.16) only at the distinct moments and only at the locations on the 
interface where the instability (tendency to form a corner) occurs. This is done as follows. 

First, we continuously check the sign of G at every marker location on the interface. 
Second, once at the certain location G becomes negative (the corner has been formed), 
we interrupt the computation, restore all variables including the local interface shape and 
curvature to their previous states at positive G, add AK 2 , A = 5/ (7o-^ 2 ) term to the nondi- 
mensional interfacial energy % v (a) at this location and solve the linear equation G(A) = 0. 
The result, 

= 15e 7 cos [4 (a + ff 7 )] - 1 

K 2 + 2Kd 2 K/da 2 + 2(dK/da) 2 ' 1 ' ' 

is deflected by 0.1 A to ensure positive G. The computation is then restarted with the regu- 
larized 7 C1) ; if needed, the procedure is repeated. This method allows to reduce the influence 
of the regularization on interface shape to a minimum by (i) applying the regularization 
selectively and (ii) choosing value for the regularization constant A no larger than needed to 
allow the computation to proceed. The maximum value A ~ 10~ 5 comes out from the com- 
putation, which translates into 5 ~ 10~ n erg; thus the additional corner energy necessary 
to regularize the problem is indeed very small. 2 

It must be noted here that the numerical method itself introduces some minor smoothing, 
since the interface is remeshed every time step with the help of cubic splines [8]. 

The regularization by curvature is not used at the contact point and its vicinity. Instead, 
to ensure stable tracking of the sensitive contact point motion over the mask edge in the 
simulations with strongly anisotropic surface energy, we had to resort to the method which 
is described in the next subsection. 

2.4.2 Contact point regularization 

e7 = 4 0) + (T- e f)tf A ( S - Sc ), (2.18) 

2 In [29], parameter 5 is introduced as the second partial derivative of the crystalline step energy with 
respect to the first coordinate derivative of step density; this is difficult, if not impossible to extract from 
experimental data. We are not aware of a single experiment there this or related quantity is measured. Since 
value is unknown, we think it makes sense to find the suitable numerical approximation to the nondimensional 
counterpart as explained. 
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where is constant value e 7 takes on the interval of the arc length < s < s c — A, T> ej^ 
is constant value e 7 takes on the interval s c + A < s < S, and H\(s — s c ) is the smoothed 
Heaviside function (see, for instance, [32]), viz. 

!0, if s — s c < —A, 
i(l + ^ + Isin^z££)), if - A<s-s c <A, (2.19) 
1, if s — s c > A. 

In (2.19), A = const, is the small half- width of the transition region from to 1 around 
s c = const. E (0, S). (2.18) and (2.19) imply that the intervals of the interface < s < 
s c — A and s c — \ < s < s c + \ evolve with lesser degree of anisotropy than the remaining 
interval s c + A < s < S. For instance, for the simulation shown in Fig. 9, we chose 
= 0.01, s c = 0.08, A = 0.02. Of course, such choice is difficult to justify physically, 
but nevertheless it allows to perform simulations and obtain the right interface shape some 
(small) distance away from the contact point [33] . 

3 Results 

Since the parameter space is very large, the overgrowth was computed for the separate cases 
of evolution with the anisotropic mobility and crystal-vapor interfacial energy. We studied 
the influence of parameters r and r y cm . Other parameters are fixed at their values cited in 
Table ??. 

3.1 Overgrowth with anisotropic mobility 

It is assumed in this section that e 7 = 0, thus ^ cv = 1 and the contact angle (f> is determined 
from the equation (2.9) where the third term identically equals zero. 

Fig. 4(a) shows the overgrowth over the blunt edge r = 0.5 (the radius = 50 nm) of the 
energetically inhomogeneous mask. Also shown in Fig. 4(b) is the actual contact angle vs. 
9. Value for 7 cm is chosen 100 (500 erg/cm 2 ), so that <p{9) > 7r/2 V 9 (partial "wetting"). One 
of accuracy checks is the comparison every time step of the actual value of the contact angle 
with the imposed (f)(9) given by (2.9). The difference is negligible for all runs performed. 
Notice that the perfectly straight facet is formed in the direction 135° to the substrate. The 
feature ("bump") is present at the junction of this facet and the horizontal, 90° facet. This 
bump is present on the surface even before the contact point reaches the mask edge, but the 
edge makes it more pronounced. The bump gradually disappears after the edge is passed. 
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Figure 4: (a) The overgrowth over the blunt edge r = 0.5 of the energetically inhomogeneous 
mask. Dashed curves show the interface at progressively increasing times, (b) The contact 
angle dependence on the angle of the normal to the mask surface, for the crystal growth 
shown in (a). 



Fig. 5(a,b) compares the overgrowth over the sharper edge r = 0.0125 (the radius = 
1.25 nm) of the energetically inhomogeneous mask as in Fig. 4(b) with the overgrowth over 
the edge of same sharpness of the energetically homogeneous mask. For the latter, value 
for the constant is taken the average of the smallest and largest values of in Fig. 4(b) 
(<f) = 1.956). The overgrowth over the energetically homogeneous mask is faster, and the 
bump is more pronounced in this case. Nevertheless, the difference in the overgrowth of the 
mask edges of different energy is quite small. 
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Figure 5: (a) The overgrowth over the sharp edge r = 0.0125; dashed curve: (f)(9) as in Fig. 
4(b), dotted curve: = 1.956. (b) The enlarged view of the tri-junction region. 



Fig. 6(a,b) compares the overgrowth over the r = 0.0125 edge of the energetically 
inhomogeneous mask that supports different (f)(9) (that is, the dynamics corresponding to 
(f)(9) > 7r/2 as in Fig. 4(b), and (f)(9) < n/2 as in Fig. 6(c)). Interestingly, the bump and 
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the crystal thickness are larger for the less wetting case (larger 0's) when the contact point 
slides over the mask edge. However, immediately after the edge the overgrowth of the more 
wetting crystal is faster, but it slows down as the overgrowth on the mask progresses. This is 
expected since (f> = n/2 on the horizontal part of the mask. Finally, the less wetting crystal 
takes over, but its thickness in the y-direction is less than one of the more wetting crystal. 



(a) (b) 
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9 

Figure 6: (a) The overgrowth over the sharp edge r = 0.0125; dashed curve: (f)(8) as in 
Fig. 4(b) (less wetting), dotted curve: (f)(9) as in part (c) of this figure (more wetting), (b) 
The enlarged view of the tri-junction region, (c) The contact angle vs. 9. To obtain this 
dependence, value for 7 cm is chosen 0.1 (0.5 erg/cm 2 ). 

The distance from the contact point to the foot point (0,0) of the mask vs. time is shown 
in Fig. 7(a,b) for the overgrowth over r = 0.5 and r = 0.0125 edge, respectively, and 
7 cm = 100. For the blunt edge, the computational data is fitted well by the cubic polynomial 
and thus the speed of the contact point vs. time is the quadratic function that has the 
minimum value Ri4x 10~ 4 (4 nm/s) at t = 125 (6.3 s) and the maximum value pa 1.4 x 10~ 3 
(14 nm/s) at t — 50,200 (2.5,10 s). For the sharp edge, the quadratic fit is better, and the 
speed of the contact point is linear function of time (the min. value ~ 3.25 x 10~ 4 (3.25 
nm/s) at t = 100 (5 s), the max. value « 5.5 x 10~ 4 (5.5 nm/s) at t = 86 (4.3 s)). 
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Figure 7: Distance vs. time (nondimensional) from the tri-junction to the foot point (0,0) 
of the mask for two values of the radius of curvature of the edge, (a) r = 0.5 (triangular 
symbols), (b) r = 0.0125 (square symbols). Fitting (solid curve) is for the time interval 
when the contact point is on the circular arc of the mask. 



This data, as well as Figures 4-6 demonstrate that de-pinning of the contact line of a crystal 
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from the solid defect of finite size does not have a threshold nature; that is, the contact line 
is in constant motion over the defect, and the problem does not have a critical eigenvalue 
such that the contact line is pinned when some parameter is smaller than this eigenvalue 
and de-pinned (moving) when the parameter is larger. In the language of [15], where the 
analogy between perturbed (on a defect) liquid contact line and an elastic spring is used 
to phenomenologically explain pinning (pages 836-38), the total pinning force exerted by 
the defect on the line is such that graphs of this force and the restoring linear elastic force 
(vs. position on a defect) are not intersecting, and thus there is no equilibrium, pinning line 
position(s) on a defect. 

Finally, the average nondimensional speed of the contact point's motion over the edge 
(naturally defined as the length of the circular segment of the mask, tcR/2, divided by the 
traversal time) is plotted in Fig. 8 vs. r. There is an evident tendency for pinning of the 
contact point at the edge as r — > 0. Data can be satisfactory fitted by a hyperbolic tangent 
curve (not shown). In the limit r = (the mathematically sharp edge) the true pinning is 
expected (the speed = 0). The contact line starts to actually "feel" the edge only when the 
latter has sufficiently small radius of curvature (r 0.1 — 0.15 10-15 nm). 

9 i 1 1 1 1 1 1 1 1 1 1 




I I i i i i i i i i i I 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



r 

Figure 8: The average speed of the contact point's motion over the mask edge. 
Note that despite such characteristic of the motion of the contact line over the edge as the 
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speed (instantaneous and average) is, of course, a function of all parameters of the system, 
the pinning on the edge as its radius of curvature approaches zero is the universal tendency. 
This was checked by means of multiple runs with different parameters. The computation 
fails for r < 0.0075 (the radius< 0.75 nm) because as the number of discretizing nodes along 
the mask increases, the nodes on the two straight segments of the mask come too close to 
each other and the interpolation becomes inaccurate. In principle, the computation with 
smaller radii can be made successful by using the nonuniform grid along the mask; we did 
not attempt to do that. It is also doubtful that very sharp mask edges (having the radii of 
curvature of the order of a molecule size (1-2 nm)) are achievable by the lithography, etching 
and polishing used in the ELO/SAG practice for the preparation of a mask. For comparison, 
the sapphire disks used for the late 1970s experiments with the liquid drops by Oliver et al. 
[16] had the edge radius of curvature smaller than 50 nm, and "much smaller" than that 
value if the defects (chipped-off hollows) appear along the edge. These disks were prepared 
by cutting from a fused sapphire boule. The aluminum disks they used had very blunted 
edges with the radius of curvature 1 to 5 /im. 



3.2 Overgrowth with strongly anisotropic ^ cv 

It is assumed in this section that a — 0, thus M = 1. The contact angle is determined from 
the equation (2.9), where the third term is nonzero due to the anisotropy in the crystal- vapor 
surface energy. 

Fig. 9 shows, as an example, the overgrowth over the blunt edge r = 0.5 (the contact 
angle <fi as in Fig. 4(b)). In Fig. 9(a), the value of the anisotropy parameter T (see (2.18)) 
is 0.07. This is just above the critical value 1/15. T = 0.085 in Fig. 9(b). For both figures, 
the value of the phase angle /3 7 is chosen 1.67 (a little larger than 7r/2). 




Figure 9: (a) The overgrowth with the strongly anisotropic interfacial energy over the blunt 
edge r = 0.5 of the energetically inhomogeneous mask. T = 0.07. (b) Same as (a), T = 0.085. 



In Fig. 9(a), the corner forms slowly at the elevation, and it is followed by the straight 
facet; other than that, the evolution resembles the isotropic evolution [8]. In Fig. 9(b), the 
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interface instantly becomes unstable and the characteristic, coarsening with time hill-and- 
valley structure is formed [29, 30]. As the result of coarsening, the single elevated corner 
(bump) followed by the a = 45° inclined facet are again formed. Thus this numerical 
experiment suggests that the strong interfacial energy anisotropy may be one of the possible 
causes for the "bumpy" crystal shapes routinely observed in ELO deposits [1, 4]. 

Fig. 10 shows the contact angle along the mask surface. In contrast to Figures 4(b) and 
6(c), this function is not monotone. As e^ -* increases quite beyond value 0.01 chosen for com- 
putations in this section, the variations in become more pronounced and the computation 
becomes unstable. 




0.2 0.4 0.6 



Figure 10: The contact angle vs. 9 for the interface evolution in Figures 9(a,b). 
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A On equation (1.1) 

The equation (1.1) follows from the following equation derived by Herring [34]: 

71 - 72 cos a 2 - 73 cos a 3 + — — sin a 2 + 77— sin a 3 = 0. (A.l) 

d(p2 d(p 3 

The equation (A.l) generalizes the Young-Dupre equation 

72 cos a 2 = 71 - 73 (A. 2) 
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to orientational dependence of surface energies (tensions) at the tri-j unction (see Fig. 11). 
To obtain (1.1) from (A.l), we identify the interface 1 in Fig. 11 with the mask-vapor 
interface, the interface 2 with the crystal-vapor interface, the interface 3 with the crystal- 
mask interface and notice that a<i = and a 3 = since the slope of the mask surface is 
everywhere continuous (Fig. 2). 




Figure 11: A sketch of the junction of three anisotropic interfaces, from [34]. 

B Parameters of the model (characteristic of epitaxy 
of GaAs-like semiconductor at 650°C ) 

This appendix has been removed from online submission to arXiv.org in order 
to comply with file size requirement; check the full version at 
www.math.buffalo.edu/~mkhenner. 
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